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Spontaneous emission and Lamb shift of atoms in absorb- 
ing dielectrics are discussed. A Green's-function approach 
is used based on the multipolar interaction Hamiltonian of 
a collection of atomic dipoles with the quantised radiation 
field. The rate of decay and level shifts are determined by 
the retarded Green's-function of the interacting electric dis- 
placement field, which is calculated from a Dyson equation 
describing multiple scattering. The positions of the atomic 
dipoles forming the dielectrics are assumed to be uncorre- 
lated and a continuum approximation is used. The associ- 
ated unphysical interactions between different atoms at the 
same location is eliminated by removing the point-interaction 
term from the free-space Green's-function (local field correc- 
tion). For the case of an atom in a purely dispersive medium 
the spontaneous emission rate is altered by the well-known 
Lorentz local-field factor. In the presence of absorption a re- 
sult different from previously suggested expressions is found 
and nearest-neighbour interactions are shown to be impor- 
tant. 



I. INTRODUCTION 

The theoretical description and experimental investi- 
gation of the interaction of light with dense atomic media 
regained considerable interest in recent years. Various 
experiments on level shifts HS], intrinsic bistability [p|0] 
and spontaneous emission gj^ in dense gases have sup- 
ported and refined the concept of local fields known for 
more than a century [[7|. Nevertheless some questions in 
this context are still not answered satisfactory even on a 
fundamental level. In the present paper I want to discuss 
one of these questions, namely the effect of an absorbing 
dielectric on spontaneous emission and level shifts of an 
embedded atom using a Green's-function approach. 

The interaction of light with dilute gases is usually 
well described in terms of macroscopic classical variables 
such as electric field and polarisation. In the macroscopic 
approach the polarisation is given by the expectation 
value of the single-atom dipole moment multiplied by 
the density of atoms [g[ . Apart from the coupling to the 
common classical radiation field, the atoms are assumed 
distinguishable and independent. This means quantum- 
statistical correlations are neglected, which is a very good 
approximation as long as the temperatures are not too 
small. It is also implicitly assumed that vacuum fluctu- 
ations of the field affect the atoms only individually and 



that the atom positions are independent of each other. 
The latter assumptions are however no longer valid in 
dense samples. 

If the resonant absorption length of some atomic tran- 
sition becomes comparable to the medium dimension d, 
i.e. for NX^d ^ 1, N being the number density and A 
the resonant wavelength, reabsorption and multiple scat- 
tering of spontaneous photons and associated effects like 
radiation trapping ||9|] or, if atomic excitation is present, 
amplified spontaneous emission need to be taken into ac- 
count. If the atomic density is further increased, such 
that iVA'^ ~ 1, one can no longer disregard the fact that 
the independent-atom approximation allows for an un- 
physical interaction of different atoms at the same posi- 
tion and Lorentz-Lorenz local field corrections are needed 
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The modification of the rate of spontaneous emission F 
by the local environment was first noted by Purcell [O . 
Alterations of this rate have been demonstrated experi- 
mentally near dielectric interfaces [[l2| , in quantum- well 
structures |13] and in cavities |14|. Based on an analysis 
of the density of radiation states Nienhuis and Alkemade 
predicted for an atom embedded in a homogeneous trans- 
parent dielectric with refractive index n [Q5|: 



F = Fon 
where Fq is the free-space decay rate. 
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p being the electric dipole moment of the transition with 
frequency cuab- The alteration of spontaneous emission by 
the index of refraction leads to interesting potential ap- 
plications as the suppression or enhancement of decay in 
photonic band-gap materials |16| . The approach of Ref. 
||l5| did neither take into account local-field corrections 
nor absorption however. 

There has been a considerable amount of theoretical 
work on local-field corrections to spontaneous emission 
of an atom in lossless homogeneous dielectrics. Essen- 
tially all approaches assume a small cavity around the 
radiating atom and the theoretical predictions depend 
substantially on the details of this local-cavity model. 
Approaches based on Lorentz's "virtual" cavity mj 18 
lead to 
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while those based on a real empty cavity [M predict 
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For pure systems or impurities in disordered non- 
absorbing dielectrics Eq.(H) is believed to be correct. On 
the other hand recent experiments with Eu'^+ ions in 
organic ligand cages verified the real-cavity expression 
Eq-CQ) (5[|§1- An explanation for the different results was 
very recently given by de Vries and Lagendijk pO| |. Ap- 
plying a rigorous microscopic scattering theory for im- 
purities in non-absorbing dielectric cubic crystals, they 
showed that the local environment determines whether 
Eq.(||) or (^ should be used. For a substitutional impu- 
rity the empty-cavity result applies, while for an intersti- 
tial impurity the virtual-cavity formula is correct. The 
latter also supports the believe that Eq. (m is the correct 
one for disordered systems like gases. 

While the effect of a transparent dielectric on spon- 
taneous emission is rather well studied, this it not the 
case for absorbing media. A first step in this direction 
was made by Barnett, Huttner and Loudon |2l|. Based 
on a discussion of the retarded Green's-function in an 
absorbing bulk dielectric they showed, that the index of 
refraction in ([^) is to be replaced by the real part n' of 
the complex refractive index n = n' + in" . They also 
argued that the square of the Lorentz-local field factor in 
(PI) should be replaced by the absolute square, leading to 



Ton'{uJab) 



n'^{uJab) + ^ 



(5) 



In order to derive this equation Barnett et al. postulated 
in p3| an operator equivalent of the Lorentz-Lorenz re- 
lation between the Maxwell and local field. This as- 
sumption has however some conceptual problems. As 
pointed out very recently by Scheel et al. p3| , an op- 
erator Lorentz-Lorenz relation cannot hold, since both 
quantities, the Maxwell field and the local field have to 
fulfil the same commutation relations. 

In a recent paper we have developed an approach that 
takes into account local-field corrections as well as mul- 
tiple scattering and reabsorption of spontaneous photons 
in modified single-atom Bloch equations [|lO[. The mod- 
ified Bloch equations provide a way of including dense- 
medium effects in a macroscopic approach. In the present 
paper expressions for the spontaneous emission rate and 
Lamb-shift of an atom in a dense absorbing dielectric or 
a gas of identical atoms are derived following the ap- 
proach of |1^. The starting point is the multipolar- 
coupling Hamiltonian in dipole approximation. The re- 
tarded Green's-function of the electric displacement field, 
which determines decay rate and Lamb shift, is calcu- 
lated from a Dyson equation in self-consistent Hartree 
approximation. As the atom positions are assumed to be 
independent from each other, local-field corrections are 
needed to remove unphysical interactions between atoms 



at zero distance. This is done in the present approach 
by an appropriate modification of the free-space Green's- 
functions rather than by introducing a cavity. The rate of 
spontaneous emission derived coincides with the virtual- 
cavity result (0) for a transparent dielectric, but differs 
from Eq.(^ in the case of absorption. It will be shown 
that in the presence of absorption near-field interactions 
with neighbouring atoms become very important, whose 
correct description requires however a fully microscopic 
approach. 



II. RADIATIVE INTERACTIONS IN DENSE 
ATOMIC MEDIA 

The present analysis is based on a description of the 
atom-field interaction in the dipole approximation using 
the multipolar Hamiltonian in the radiation gauge [p4| 
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Here dj is the dipole operator of an atom at position 
rj. D is the operator of the electric displacement with 

V • Z? = 0. 

It was shown in jlO) that the effects of radiative atom- 
atom interactions in a dense medium can be described in 
Markov approximation with a nonlinear density-matrix 
equation 



Hq,p 



+ i 



Pi^ 



'^t^^Ln+'^i^tn-^P 



'^'^LiV 



<^l(^tJ,,P 



ih'r 



<,crp 



-I^^^t 



"2~ 



\ <^l<^f^P + P<^l<^f^ - 2cr^pcrJ I 



c^P' [^i'P] + '^i' [f^P'P] 



(7) 



Here p^ is the dipole matrix element for a polarisation 
direction e^ and cr^ , aj^ are the corresponding atomic low- 
ering and razing operators. The first term describes the 
free atomic evolution and the second the interaction with 
some local classical field £l- hfj,i, and F^;^ are matrices, 
whose eigenvalues yield Lamb shifts of excited states and 
spontaneous emission rates. V^ and /i^^ describe collec- 
tive relaxation rates and light-shifts due to the incoherent 
background radiation generated by absorption and ree- 
mission of spontaneous photons (radiation trapping) . 

It should be noted that the incoherent background ra- 
diation causes a decay as well as an incoherent excitation 
with equal rate F'^. Thus F"^, which is proportional to 
the excitation of the host medium |10(| , describes induced 
mixing processes, while F can be interpreted as the rate 
of spontaneous decay. Similarly h'^ describes a light-shift. 



which for a two level system is equal in strength and op- 
posite in sign for the ground and excited state. It is also 
proportional to the excitation of the host medium and 
can thus be interpreted as induced light shift. In con- 
trast ft, is a frequency shift of an excited state only and 
does not require excitation of the host medium. 
The matrices T^i, and h^i, are given by pO[| 
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where D^^(af, cj) = J_^ dT_D^i,(a;, r)e'"'^ is the Fourier- 
transform of the retarded Green's-function (GF) of the 
electric displacement field defined here as 



D^,{x,t)=9{t){0\ D^,{ruh),D,{r2,h) |0) 



(10) 



with X — fi — f2 and t — ti^t2. In the case of randomly 
oriented two-level atoms, one can replace p^ —^ p and 
perform an orientation average yielding a single decay 
rate F and a single excited-state level shift h. 

The dense atomic medium affects the spontaneous 
emission of a single probe atom due to multiple scattering 
of virtual photons. The scattering process can formally 
be described by a Dyson equation for the exact retarded 
GF 



The over-bar denotes an average over some possible in- 
homogeneous distribution and N is the number density 
of atoms. 

With the above made approximations, the Dyson 
equation (nw contains also scattering processes between 
atoms at the same position. In a continuum approxima- 
tion the probability of two point dipoles being at the same 
position is of measure zero. This nevertheless leads to a 
non- vanishing contribution, since the dipole-dipole inter- 
action has a (5-type point interaction. This unphysical 
contribution needs to be removed by a local-field correc- 
tions, which will be discussed in the following section. 



III. LOCAL-FIELD CORRECTION OF 

FREE-SPACE GREEN'S-FUNCTION AND 

LORENTZ-LORENZ RELATION 

The retarded Green's-function in 

free space D0,(1, 2) = e{h-t2){0\ [D°^{1) , D°{2)] \0)eo^ 
where 1; 2; . . . stand for ri, ii; r2, h', ■ ■ ■ etc., is a solution 
of the homogeneous Maxwell equation with i5-like source 
term 
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D(l, 2) ^ D0(1, 2) - // .3 .4 D0(1, 3) n(3, 4) D(4, 2). D" has a particularly simple form in reciprocal space ^ 



(11) 

Here the integration is over t from — oo to -l-oo and the 
whole sample volume. D'' is the (dyadic) GF in free space 
and n is a formal (dyadic) self-energy. As shown in [Q , 
the self-energy can be described for randomly oriented 
two-level atoms in self-consistent Hartree approximation 

by 



n(i=2)^^|f^0(ii-t2)(H(ii),a,(t2)] 
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where k = lo/c, and A^ = 1 — ^^. It may be worthwhile 
noting, that D'' is not transverse in g-space, although Vi • 
D°(l,2) = 0. The corresponding function in coordinate 
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1 is unity matrix and <t = |6)(a| is the atomic spin- flip 
operator from the excited state \a) to the lower state \b) 
in the Heisenberg picture, i.e. it contains all interactions. 
The factor 2/3 results from an orientation average. 

We now make a continuum approximation and assume 
a homogeneous medium, such that 
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Here x = lafl and 



P{z) = l-- + ^, 



Q(z) = -l + --4. (19) 



where 



n(l,2)^p(ii,t2)(S(fi-r2)l, 



(13) 



P(ii,t2) = 3|2^^(ii-*2)([at(^i),c7(t2)]). (14) 



One recognises from Eq.(|lq) that the retarded GF of 
the dipole-dipole interaction contains a 5-type point con- 
tribution. In order to eliminate the unphysical interac- 
tions between different atoms at the same position, one 
has to remove this term from the GFs in the scattering 
part of the Dyson equation (|lj) . 



D°(f,cj+) -^ F"(f,cj+) = D"(.f,cj+) - ^—5{x)l. (20) 

With this local-field corrections we obtain a modified 
Dyson equation (in reciprocal space) 
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and introducing F{q,Lo) = D((f, w) — ih/Seo 1 we arrive 
at 

F{q,oj+) = F°{q,oj+) - F"{q,cu+)p{uj+)F{q,uj+). (22) 

In reciprocal space one finds 



F°(g,c.+) = D°(<r,c.+)-^l, 
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Eq.(b2|) can easily be solved to yield 
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where we have introduced the dynamic polarisability of 
the atoms 
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The poles ±go of Eq.(E5p determine the (in general non- 
linear) complex dielectric function 
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This is the well-known Lorentz-Lorenz relation between 
the microscopic polarisability a and the complex dielec- 
tric function e{oj). Thus we have shown that the local- 
field correction of the free-space Green's-function ( |20| ) is 
exactly the one that reproduces the well-known Lorentz- 
Lorenz relation. 



IV. MODIFICATION OF SPONTANEOUS 
EMISSION AND LAMB-SHIFT 

Eq. (E5|) can be transformed back into coordinate space 
using F{x,uj+) = (27r)-3/d3gF(g,a;+)e-^«~. The 
Fourier-transform of the projector {qoq) yields spher- 
ical Bessel functions [g5[. For the present purpose we 
however need only the orientation-averaged quantity 
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One recognises that the Fourier-transform of F((f, u;+) 
diverges for x ^ 0, which is due to the large-g behaviour 
of the GF. In order to remove these singularities one can 
modify the GF by introducing a regularisation. Physi- 
cally the singular behaviour at a; — *■ is due to the fact 
that atoms very close to the atom under consideration 
can have a large effect on spontaneous emission and level 
shifts. One cannot expect the continuum approxima- 
tion used here to yield accurate results on length scales 
comparable to the mean atom distance. Here rather a 
fully microscopic description of very close atoms includ- 
ing their motion (collisions) is needed. This is however 
beyond the scope of the present paper and we therefore 
restrict the analysis to a regularisation of the Green's- 
function. There is no unique regularisation procedure, 
and we here just choose a convenient one 
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With this we find in the limit A » |go| 
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where R = fc/(\/2A). It is important to note, that F is 
exactly causal, iie{oj) fulfils the Krames-Kronig relations. 
This would not have been the case if as according to 
the result of Barnett et al. |2^J22] the absolute square 
|(£ + 2)/3|2 would be present instead of ((e + 2)/3)2. 

With this result we find for the decay rate and excited 
state Lamb-shift 
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For an atom in a purely dispersive disordered medium, 
i.e. for e" = 0, the second term in Eq.(|3l|) for the spon- 
taneous decay rate vanishes identically and we are left 



with the "virtual" cavity result Eq.(g). Likewise there 
are no contributions from the first term in ( |32|) to the 
Lamb shift in this case. 

In the presence of absorption, that is if the probe-atom 
transition frequency comes closer to a resonance of the 
surrounding material (as it would naturally be the case 
for a collection of identical atoms) F is different from 
the result obtained in ||2l|,^.g8[. In this case there are 
also non-vanishing terms that contain the regularisation 
parameter R~^ and R^^. These terms must be inter- 
preted as contributions due to resonant energy transfer 
with nearest neighbours, a process which cannot accu- 
rately be described in the continuum approach used here. 

As the Lamb shift is concerned, Eq.(p2|) shows that in 
a purely dispersive medium, that is far away from any 
resonances only nearest-neighbour interactions matter. 
This is intuitively clear since in this case the transition 
frequency is only affected by collisions. Only in the pres- 
ence of absorption there is also a bulk-contribution to the 
Lamb shift as described by the first term in (p2[). 

For a dense gas of identical atoms or of atoms of the 
same kind but with some inhomogeneous broadening, 
Eqs.(pi|) and ( [3^ ) are only implicit, since the complex 
polarisability e depends on the decay rate and level shift. 
Hence a self-consistent determination of F and h is nec- 
essary. If the density of atoms is much less than one 
per cubic wavelength one can consider an expansion of 
F and h in powers of the atomic density N. Defining 
a = a' + ia" one finds with Eq. ( p7| ) for the bulk contri- 
butions 
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In the case of radiatively broadened two-level atoms, the 
real part of atomic polarisability vanishes at resonance, 
i.e. a' — 0. Thus in lowest order of the density there is 
only a contribution to the excited state frequency pro- 
portional to the population difference between excited 
and ground state. For an inverted population the tran- 
sition frequency is red-shifted, for balanced population 
the level shift vanishes and for more atoms in the lower 
state the transition frequency is blue shifted. As a re- 
sult spontaneously emitted radiation from an initially in- 
verted system will have a chirp very similar to the chirp 
in Dicke-superradiance |g6[ . It should also be mentioned 
that the shift of the transition frequency discussed here 
is physically different from the familiar Lorentz-Lorenz 
shift. The LL-shift is due to the dispersion of the in- 
dex of refraction at an atomic resonance and is thus in 
contrast to the absorption a" independent on Doppler- 
broadening pTI]. 



V. SUMMARY 

In the present paper we have discussed the rate of spon- 
taneous emission and the excited-state level shift of a 
two-level type probe atom inside a homogeneous, disor- 
dered absorbing dielectric. The dielectric was modelled 
by a collection of atomic point dipoles, which also in- 
cludes the case of a dense gas of identical atoms. The 
multiple scattering of photons between the atoms (dipole- 
dipole interaction) was described by a Dyson integral 
equation for the exact retarded Green's-function of the 
electric displacement field in self-consistent Hartree ap- 
proximation. The atoms were assumed distinguishable 
with random independent positions. The latter assump- 
tion made a continuum approximation possible and the 
Dyson equation could be solved analytically. In order 
to exclude unphysical dipole-dipole interactions of dif- 
ferent atoms at the same position arising in the con- 
tinuum approximation with independent atomic posi- 
tions, a local-field correction of the free-space retarded 
Green's-function was introduced. This lead to the well- 
known Lorentz-Lorenz relation between the complex di- 
electric function e(w) and the nonlinear atomic polaris- 
ability a{uj). The expression for the spontaneous-decay 
rate found by this method agrees with the virtual cav- 
ity result ||l^ in the absence of absorption. This is 
an expected result for atoms in disordered dielectrics 
pO{ . It was shown that the excited-state Lamb shift is 
in this case only affected by nearest-neighbour interac- 
tions, which could not be treated accurately within the 
present approach however. In the presence of absorption 
the spontaneous-emission rate differs from the results ob- 
tained in |pi| , p2y28| in two ways. First there are impor- 
tant nearest-neighbour contributions, which were absent 
in the models of pi| , p2|. Secondly the bulk-contribution 
is different form Refs. ]2l| , |2^ , p8| | , since causality of the ex- 
act retarded GF requires the Lorentz-field factor to enter 
as square and not as absolute square. It is interesting to 
note, that apart from a small difference in the R~^ term, 
the decay rate derived here is identical to one very re- 
cently obtained by Scheel and Welsch ||2^] on the basis of 
a completely different approach, namely a quantisation 
of the electromagnetic field in a linear dielectric. 



ACKNOWLEDGEMENT 

I would like to thank Charles Bowden, Janne Ru- 
ostekoski and Dirk-Gunnar Welsch for stimulating dis- 
cussions and D.G. Welsch and S. Scheel for making Ref. 
p9| available prior to publication. 



[1] V. A. Sautenkov, H. van Kampen, E. R. Eliel, and 

J. P. Woerdman, Phys. Rev. Lett. 77, 3327 (1996). 
[2] H. van Kampen, V. A. Sautenkov, C. J. C. Smeets, E. R. 

Eliel and J. P. Woerdmann, Phys. Rev A 59, 271 (1999). 
[3] J. J. Maki, M. S. Malcuit, J. E. Sipe and R. W. Boyd, 

Phys. Rev. Lett. 67, 972 (1991). 
[4] M. P. Hehlen, H. U. Gudel, Q. Shu, J. Rai, S. Rai, and 

S. C. Rand, Phys. Rev. Lett. 73, 1103 (1994). 
[5] G. L. J. A. Rikken, and Y. A. R. R. Kessener, Phys. Rev. 

Lett. 74, 880 (1995). 
[6] Frank J. P. Schuurmans, D. T. H. de Lang, G. H. Weg- 

dam, R. Spirk, and A. Lagendijk, Phys. Rev. Lett. 80, 

5077 (1998). 
[7] H. A. Lorentz, Wiedem. Ann. 9, 641 (1880); L. Lorenz, 

Wiedem. Ann. 11, 70 (1881); L. Onsager, J. Am. Chem. 

Soc. 58, 1486 (1936); C. J. F. Bottcher, Theory of electric 

polarization (Elsevier, Amsterdam, 1973); M. Born and 

E. Wolf Principles of Optics, (Wiley, New York, 1975); 

J. van Kronendonk and J. E. Sipe, in Progress in Op- 
tics, ed. by E. Wolf (North-Holland, Amsterdam, 1977) 

Vol. XV; J. T. Manassah, Phys. Rep. 101, 359 (1983); 

C. M. Bowden and J. Bowling, Phys. Rev. A 47, 1247 

(1993); ibtd 49, 1514 (1994); O. Morice, Y. Castin, and 

J. Dalibard, Phys. Rev. A 51, 3896 (1995); J. Guo, A. 

Gallagher, and J. Cooper, Opt. Comm. 131, 219 (1996); 

J. Ruostekoski and J. Javanainen, Phys. Rev. A 56, 2056 

(1997), ibid. 55, 513 (1997). 
[8] see for example: M.Sargent III, M.O. Scully, and W.E. 

Lamb, Laser Physics (Addison Wesley, Reading, MA 

1974). 
[9] T. Holstein, Phys. Rev. 72, 1212 (1947); ibid. 83, 1159 

(1951); 
[10] M. Fleischhauer and S. F. Yelin, Phys. Rev. A m press 

(March 1999). 
[11] E. M. Purcell, Phys. Rev. 69, 681 (1946). 
[12] K. H. Drexhage, in Progress in Optics, ed. by E. Wolf 

(North-HoUand, Amsterdam 1974), Vol. XIL 
[13] Y. Yamamoto, Opt. Comm. 80, 337 (1991). 
[14] F. De Martini et al. J. Opt. Soc. Am. B 10, 360 (1993). 
[15] G. Nienhuis and C. Th. J. Alkemade, Physica 81C, 181 

(1976). 
[16] E. Yablonovitch, Phys. Rev. Lett. 58, 2059 (1987); E. 

Yablonovitch and T. J. Gmitter, ibid 63, 1950 (1989); S. 

John and T. Quang, Phys. Rev. A 50, 1764 (1994). 
[17] J. Knoester and S. Mukamel, Phys. Rev. A 40, 7065 

(1989). 
[18] P. W. Milonni, J. Mod. Optics 42, 1991 (1995). 
[19] R. J. Glauber and M. Lewenstein, Phys. Rev. A 43, 467 

(1991). 

P. de Vries and A. Lagendijk, Phys. Rev. Lett. 81, 1381 

(1998). 

S. M. Barnett, B. Huttner, and R. Loudon, Phys. Rev. 

Lett. 68, 3698 (1992). 

S. M. Barnett, B. Huttner, R. Loudon, and R. Matloob, 

J. Phys. B 29, 3763 (1996). 

S. Scheel, L. Knoll, D -G. Welsch, and S . M. Barnett, 

(unpublished; preprint |quant-ph/9811067[ 24.Nov.1998). 

see for example: J. R. Ackerhalt, and P. W. Milonni, J. 



[2o; 

[21 

[22; 
[23; 
[24; 
[25; 



Mod. Phys. 70, 447 (1998). 
[26] see for example the review on superradiance: M. Gross 

and S. Haroche, Phys. Rep. 93, 301 (1982). 
[27] J. Guo, A. Gallagher, and J. Cooper, Opt. Comm. 131, 

219 (1996). 
[28] G. Juzelitinas, Phys. Rev. A 55, R4015 (1997). 
[29] S. Scheel und D. G. Welsch, unpublished 



Opt. Soc. Am. B 1, 116 (1984). 

P. de Vries, D. V. van Coevorden, and A. Lagendijk, Rew. 



